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Abstract. A detection of primordial non-Gaussianity could transform our understanding of 
the fundamental theory of inflation. The precision promised by upcoming CMB and large- 
scale structure surveys raises a natural question: if a detection given a particular template is 
made, what does this truly tell us about the underlying theory? 

In this paper we present a systematic way to constrain a wide range of non-Gaussian 
shapes, including general single and multi-field models and models with excited initial states. 
We present a separable, divergent basis able to recreate many shapes in the literature to 
high accuracy with between three and seven basis functions. The basis allows shapes to be 
grouped into broad "template classes", satisfying theoretically-relevant priors on their diver- 
gence properties in the squeezed limit. 

We forecast how well a Planck-like CMB survey could not only detect a general non- 
Gaussian signal but discern more about its shape, using existing templates and new ones we 
propose. This approach offers an opportunity to tie together minimal theoretical priors with 
observational constraints on the shape in general, and in the squeezed limit, to gain a deeper 
insight into what drove inflation. 
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1 Introduction 

An early period of accelerated expansion, perhaps a trillionth of a second after the Big Bang, 
is proposed to solve a number of problems unresolved by the Big Bang scenario, such as the 
flatness and horizon problems. The paradigm of single- field slow-roll inflation is the simplest 
model to describe this acceleration, and makes broad predictions of adiabatic and Gaussian- 
distributed primordial density (scalar) perturbations, described by a nearly scale-invariant 
2-point correlation, and smaller gravitational metric (tensor) perturbations. In this model, 
the scalar and tensor amplitudes and scale dependence are related through a 'consistency 
relationship'. 

The precision of astrophysical measurements has dramatically improved over the last 
decade. Cosmic Microwave Background (CMB) measurements, such as from the Wilkin- 
son Microwave Anisotropy Probe [1, 2], small-scale CMB measurements from the Atacama 
Cosmology Telescope and South Pole Telescope [3, 4], and large-scale structure (LSS) ob- 
servations such as from the Sloan Digital Sky Survey [5, 6], are all entirely consistent with 
single- field slow-roll predictions, placing strong constraints on the scalar power spectrum and 
upper limits on the degree of deviations from Gaussianity and amplitude of tensor modes. 

The agreement between single-field inflationary predictions and observations is a pro- 
found success for cosmology, but it is as yet, insufficient to inform us about the underlying 
theory from which inflation derives. Rapid theoretical progress in high energy effective field 
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theory has led to a wide range of possible Lagrangians for inflation [7-9]. These often go be- 
yond making distinct predictions for the form of the single- field inflationary potential and can 
include multiple dynamical fields and non-canonical derivative (kinetic) terms in the action. 

These alternative mechanisms can generate new observational signatures, including dif- 
ferent consistency relationships relating scalar and tensor perturbations [10, 11], the addition 
of non-adiabatic (isocurvature) modes [12, 13], and the possibility of observationally mea- 
surable non-Gaussian correlations [14]. In particular, the possibility that primordial non- 
Gaussianity may be detectable as a non-zero 3-point correlation function, or bispectrum, has 
been a major development in the search for observational signatures of the underlying infla- 
tionary theory. What is most exciting is that different theories can give rise to bispectra with 
distinct scale dependencies, such that measuring not only the amplitude but also the scale- 
dependence, or 'shape', of the bispectrum could provide a direct insight into the inflationary 
mechanism. 

Much work has focused on potentially measuring the amplitude of commonly predicted 
shapes, such as the local [14-16], equilateral [17], and orthogonal [18] templates. Recent 
theoretical developments have also led to a wider population of bispectra, including those 
from fast-roll inflation [19-22], quasi-single field inflation [23, 24], warm inflation [25, 26], and 
excited initial states [27-29]. There are also hybrids of multi- field and non-slow-roll models 
[30-32], and the inclusion of isocurvature modes in the non-Gaussian correlations [33-35]. 
These bispectra can have very different shapes, meaning their signal is weighted towards 
different configurations of the 3 wavenumbers in (Fourier) /c-space. How divergent shapes 
are in the 'squeezed' /c-configuration, when one of the three length scales contributing to the 
3-point function becomes much larger than the other two, in particular can signal whether 
inflation is derived from a single-field or multi-field model. 

The divergence in the squeezed limit could also be constrained by its effect on large-scale 
structure. A non-Gaussian signal peaking in the squeezed limit would directly couple large 
scale modes to small scales, on which non- linear halos are forming [36]. This gives rise to 
an additional contribution to the halo bias, determining how the number density of halos 
of a given mass are related to the underlying linear power spectrum. In theory, wide field 
large-scale structure surveys could provide a sensitive constraint on the divergence properties 
of non-Gaussianity [37-43]. 

Given the diversity of theoretically motivated shapes, an intriguing question is how well 
one might actually be able to determine the shape of primordial non-Gaussianity, rather than 
purely assuming a shape template is the true shape a priori. To what extent can the shape of 
non-Gaussianity be reconstructed using the CMB and LSS 3-point correlations? If a positive 
detection is made assuming a template, how well would such a detection really constrain the 
underlying shape and the theoretical model that generated it? 

This 'reconstruction' approach has been widely considered in the context of the infla- 
tionary power spectrum, both in terms of P{k) reconstruction (e.g. [44-46]), and measuring 
the hierarchy of slow-roll parameters (e.g. [47-50] and references therein), instead of assuming 
a nearly scale-invariant spectrum parametrized by a constant tilt and a constant running 
dus/dlnk. 

Unfortunately, calculating theoretical predictions for CMB bispectra is computationally 
cumbersome in its exact form, requiring 4-dimensional integrals to be performed. A formalism 
to make the calculation tractable for general bispectra was introduced in [51]. The authors 
proposed a technique to create templates for shapes by expanding non-separable shapes on 
a basis set of bispectra that are explicitly separable functions of the three wavenumbers. 
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The separability reduces the 4-dimensional integral to a tractable computation without a 
significant reduction in the accuracy of the computed CMB bispectra. This approach has 
been used to forecast constraints for a variety of fundamental shapes [2] and adapted to other 
basis sets to describe oscillatory, rather than monotonic, shapes [52]. 

In this work we present an alternative separable basis to efficiently describe and investi- 
gate the broad class of nearly scale-invariant general bispectra in terms of their squeezed limit 
properties. We discuss a way to expand a general shape in the basis, which is specifically 
tuned to enable us to systematically increase the complexity of the template in a theoretically 
motivated way. We forecast the potential for determining the underlying non-Gaussian shape 
given upcoming CMB temperature and E-mode polarization data modeled on the Planck 
survey. 

The format of the paper is as follows. In section 2, we review the formalism used to 
calculate CMB bispectra. We introduce a separable basis to describe general shapes that 
are scale-invariant and potentially divergent, and discuss how this basis can be applied to 
describe a wide variety of shapes in the literature. Using the basis, we develop an expansion 
that allows us to incrementally investigate classes of bispectra motivated by theories. In 
section 3, we present a Fisher analysis quantifying how well a Planck-like survey will be 
able to distinguish between and constrain individual bispectrum shapes. Using a principal 
component analysis, we find the best to worst measured uncorrelated shapes, and compute 
the overall uncertainties in the bispectrum measurement as a function of /c-space configuration 
under different theoretical priors. We use these results to establish how much we can learn 
about the bispectrum shape, and hence with what confidence we might be able to narrow 
down the underlying inflationary theory. In section 4, we summarize our findings and discuss 
implications for future work. 

2 Efficient calculation of a general non-Gaussian shape 

In this section we lay out the formalism to describe and compute general bispectra. In sub- 
sections 2.1 and 2.2, we respectively review the calculation of the CMB bispectrum given the 
primordial 3-point function and the definitions of covariances in wavenumber and multipole 
space that roughly quantify the theoretical similarity of two bispectra. In subsection 2.3 we 
introduce a separable basis set to describe general bispectra and develop computationally 
tractable templates. Subsection 2.4 discusses the application of the basis set to a variety of 
theoretical bispectra and templates in the literature. How bispectra can be classified and 
presented pictorially is reviewed in subsection 2.5. 

2.1 The CMB bispectrum 

While Gaussian fluctuations are wholly described by a 2-point correlation function, a full 
description of non-Gaussian fluctuations requires higher order correlations that are not triv- 
ially related to the 2-point function. The simplest higher order correlation is the 3-point 
function, where the 3-point Fourier space statistic analogous to the 2-point power spectrum 
is the bispectrum, B^, defined by 

($(ki)$(k2)«'(k3)) = {27rfS{k^ + k2 + k3)S$(fcl, k2, h). (2.1) 

$(k) is the primordial gravitational potential, related to the curvature perturbation by 
$ = |7^. Under the assumptions of statistical isotropy and homogeneity, the bispectrum 
is dependent on only the magnitudes of the wavenumbers, ki, k2, and k^. 
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The bispectrum is often parameterized by a shape, ^(A;!, ^2, k^), and an amphtude, /nl, 
at an arbitrary configuration in /c-space which together determine the bispectrum at all scales, 

(''^''^''^)^_B^{h,k2,k3) = fNLS{ki,k2,ks). (2.2) 



The typical convention is to choose A'' = 6[27r^ (|)^ A^(A;o)]^, where A^(fco) is the amplitude 
of the primordial power spectrum of the curvature perturbations at a pivot scale /cq. Shapes 
are typically normalized such that S{ko, kQ, /cq) = 1- 

CMB statistics are commonly described by correlations between angular moments on 
the sky, aimi calculated through a spherical harmonic decomposition of the photon transfer 
functions, Ai{k), integrated along the line of sight and sourced by the primordial perturba- 
tions, 

j^A,{kMk)Y,M. (2.3) 
The CMB 3-point correlation function is given by 



2X3 



{ai^miae^miaeams) = ( " ) / dxx'^ / dkidk2dkz{kik2k-ifBq,{ki,k2,kz) 



X (^i)^^2 (^2) {k^)je^ {kix)ji2 {k2x)je., (k^x) 

X j dOxY^^mi(x)Y£2 m2(^)-^3m3(^)- (2-4) 

To perform the integrals over k and x, we use the CAMB^ code [53], which uses the line of 
sight approximation [54] to calculate the photon transfer functions Ag. 

We consider purely isotropic bispectra for which the integral over J df^x is a separable 
geometrical factor called the Gaunt integral. The properties of the Gaunt integral require 
that the non-zero correlations have £1, £2, and ^3 satisfying an even sum ii + £2 + £3 and 
\£i — £21 < £3 < £1 + £2 for £1, £2 < £3- Under the assumption of isotropy, the angle- aver aged 
angular bispectrum is the 3-point analogue to the Ci, 

Bhiih = X] ( rl\ m2 ma ) ^"^i™i"^2m2«^3'n3), (2-5) 

where the bracketed term is the Wigner-3j symbol. To further separate out a purely geomet- 
rical factor from the angular- aver aged bispectrum, it is convenient to work with the reduced 
bispectrum hg^i^i^^ 



so that 



/ (2^1 + 1) (2^2 + 1) (2^3 + 1) (£x£2h \, 

^^ife£3 - V I ' ^^^2^3' ^ 



— ^ y dxx'^ J dkidk2dk3{kik2k3)'^B^{ki,k2,k3) 
X A^^ {ki)Ai^ (^2) {k3)ji^ {kix)ji^ {k2x)ji>., (ksx). (2.7) 



http://camb.info 
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Here x is a dummy variable which should be integrated between zero and infinity. We note that 
in an analogous evaluation of C^, x has a physical interpretation as the comoving distance 
to the surface of last scattering. One might assume therefore that in the 3-point integral 
the upper limit Xmax could be set to (tq — Tree)- However, as others have previously also 
commented [18], the integral over x arises out of rewriting the delta function in (2.1) as an 
integral over a product of Bessel functions. Numerically, for a general bispectrum we find the 
value of Xmax ensuring the required degree of convergence is ^-dependent, and typically needs 
to be greater than (tq — Tr, 



' recj ■ 



2.2 Shape similarity 

The degree to which a bispectrum B is theoretically similar to another, B' , can be quantified 
by a A;-space correlation coefficient, or 'A;-space cosine', corr^, integrated over the fc-space 
tetrapyd volume V with weight w [51, 55], 

{S,S') = / S{ki,k2,ks)S'{ki,k2,k3)w{ki,k2,ks)dkidk2dk3, (2.8) 



corrk(5,S') = ^iM2^=. (2.9) 
^{S',S')k{S,S)k 

An analogous statistic describing the similarity of two bispectra in multipole space can be 
quantified by their £-space correlation coefficient, or '^-space cosine', corr^ [51], 

corre{B,B') = ^JM2^^. (2.11) 
^{S,S)e{S',S>), 

These correlation statistics are frequently used within non-Gaussian shape studies to 
quantify how well a template matches a given shape. A template can be obtained, given a 
set of n basis shapes by applying the Gram-Schmidt algorithm to give a set of orthonormal 
basis functions {TZn} (in either k or i space) that can be used to create a matched template 
for a specific non-separable shape S, 



S^Xlate = E^^- ^)^^ = E (2.12) 

i=0 



The efficacy of the template can then be quantified by the cumulative cosine, corr(S', 'S'l^^p/ate) 
as in (2.9) or (2.11). 

A high correlation coefficient signals a good fit. If this cosine is close to one then the 
two shapes are sufficiently alike that one might expect constraints on the amplitude of B 
can be taken as constraints on i?' as well, without having to do a separate analysis of the 
data. If on the other hand, the cosine is low, then it is likely that separate analyses of the 
data are needed for B and B\ because a template for B will not be able to pick out a non- 
zero signal for B\ and vice versa. We note though that while a template may have a large 
cosine with the shape, this does not automatically mean that the template will be able to 
accurately model a correlation between the true shape and a third shape, with which it is not 
similar. The extreme example of this would be if the third shape were exactly proportional 
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to the discrepancy between the template and shape. In constructing a template, if this was a 
concern, one might want to tailor it to the purpose by altering the weight in the Gram-Schmidt 
decomposition to ensure a minimization of the covariance between the shape and template 
over a given region of (k or i) configuration space in which a third shape was relevant. 

While it is extremely useful to establish the similarity of shapes, it is the converse of 
this, how well two shapes can be distinguished from one another, using data, that is the 
main focus of this work. This provides a motivation to consider an efficient way to generate 
-£-space bispectra explicitly by creating templates described by basis functions separable in 
ki, k2 and ks as we discuss below. To do this, in sections 2.3 and 2.4 we develop a framework 
to describe the possible degrees of freedom that a general shape might have under a specific 
theoretical prior. Furthermore, neither corr^ nor corr^ takes into account the noise, sky cover- 
age or resolution characteristics of a particular survey. A Fisher or analysis that includes 
the measurement covariance matrix essentially uses a refined, survey-dependent extension of 
(2.11) that tailors the correlation statistic to reflect the observational, rather than intrinsic, 
distinguishability of shapes. In section 3, we perform a Fisher analysis and use the correlation 
statistics to quantify how well upcoming surveys might distinguish one shape from another. 

2.3 A new separable basis, /C„ 

The 4-dimensional integral over the product of highly oscillatory functions given in (2.7) is 
computationally intensive. This has been a barrier to efficiently calculating observational 
predictions for the CMB bispectrum. As a result many studies have focused on models for 
which the primordial bispectrum can be written as, or well-approximated by, a separable 
(symmetric) function of (ki,k2,k^), 

S{ki,k2,k3) = fiki)g{k2)h{k3) + cyclic perms, (2.13) 

such that the 3-dimensional integral over ki,k2, and in (2.7) is reduced to a product of 
three 1-dimensional integrals. 

In [51] the authors proposed a way to reduce the computation time for general models 
by expanding the shape in terms of a separable basis, 

S{kl, k2, k3) = J2 (^^Qi{kl,k2, k3). (2.14) 

i 

Each Qi is constructed from symmetrized products of three 1-dimensional polynomials of ki, 
k2, and k^, 

p r s 

Qniki, k2, kz) oc ^ Cpi ^ Crj ^ Csk (^klk^k^ + k{k^kl + 4 perms) , (2.15) 

i=0 j=0 k=0 

where n maps onto a combination of {p,r,s} > 0, and Cpi, Crj, Cgk are constants. 

Using the Q„ basis, the equilateral template can be reconstructed to 98% accuracy (ac- 
cording to the cumulative cosine) using 6 basis functions [51], while other shapes motivated by 
single-field inflation models can require 20 or more mode functions to get > 95% convergence 
[56]. 

An analysis of data to constrain the bispectrum depends not only on the uncertainties 
inherent in the data itself, but also the theoretical priors determining the model being com- 
pared with the data. The choice of a separable basis set to describe the theory is therefore 
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also influenced by this prior. An analysis allowing the primordial bispectrum to take any form 
(i.e. with no shape prior on the forms of the separable functions /j, gi, hi) would use a discrete 
set of fc-space bins to describe the uncorrelated amplitudes at each scale and configuration. 
In such a scenario, no theoretical prior is applied and the constraints on the bispectrum are 
simply those determined by the data. In studying theoretically motivated models of infla- 
tion, however, there can be broad or specific characteristics of the bispectra that it would 
be reasonable to impose in conjunction with the data that suggest a form for the separable 
basis functions. For example, a Fourier basis may be more efficient than a polynomial one for 
describing bispectra with oscillatory features [52]. The two minimal assumptions we consider 
here as theoretical priors are that the bispectrum i) has a roughly monotonically changing 
amplitude as a function of scale, and ii) like the power spectrum, it is nearly scale-invariant. 

The polynomial basis of [51] does not naturally confine shapes to these two common 
theoretical properties of bispectra. Firstly, the polynomial basis does not naturally restrict 
itself to scale-invariant shapes, because i + j + k > in (2.15); resulting sums of the basis 
functions are thus scale-dependent in general. Most theoretically-motivated bispectra in the 
literature, however, are nearly scale-invariant, with i + j + k ^ (see [57] for a review). Such 
shapes can be reduced to functions of two variables, ki/k^ and ^2/^3- There are exceptions 
to this, of course, such as non-Gaussianity from particle production [58] or from features in 
the inflationary potential [59, 60]. These models can strongly deviate from sc ale- invar iance 
because modes leaving the horizon at a specific moment, when particle production is occurring 
or a feature in the potential is important, are preferentially populated. Secondly, different 
types of theoretical mechanisms generating bispectra predict different divergence properties 
in the squeezed limit, where fei <C A:2 ~ ^3. We consider the squeezed limit as ki = ke, 
the long wavelength mode, and k2 = k^ = kg, the short wavelength modes, so that for 
scale-invariant shapes the squeezed limits purely dependent on Xsq = ki/kg. Single- field 
inflation models, through a consistency relation [61] predict the bispectrum will vanish in this 
limit. Local bispectra, typically arising in multi-field models, have a x'^^ divergence, while 
excited states can have divergence. Since the powers i, j, and k in (2.15) are > 0, the 
{Qn} all tend towards a constant value in the squeezed limit, and thus cannot effectively 
describe shapes diverging in the squeezed limit. As a result, this basis is not immediately 
suited to reconstructing templates that display specific divergence behaviors in this limit, 
without a further prior being imposed. For example, the compelling and well-studied local 
template cannot easily be recovered using the {Qn} basis without either using many more 
basis functions or ignoring the divergent part of the shape that makes the local template 
distinct from others [56]. 

In this paper we introduce a set of separable basis functions, {/Cn}, that efficiently 
describe nearly scale-invariant and potentially divergent shapes, and explicitly consider the 
forms of the shapes generated using this basis under various divergence constraints. Explicitly, 
we consider 



Here f^j^ are expansion coefficients and n again denotes a combination of powers {p, r, s} of 
the wavenumbers {ki,k2,ks). /C„ is defined as 




(2.16) 



n 



ICn{ki,k2,k3) 



1 



^1 ^2 ^3 + {prs} perms 



(2.17) 
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where Mn is the number of distinct permutations of {p, r, s}. p' is defined as 

p=2+^-^ )\ iZ^ (2.18) 

and similarly for r' and s' . 

Equations (2.17)-(2.18) ensure that each /C„ is normalized in the conventional way, with 
ICn{ki,k2,k3) = 1 at ki = k2 = = ko. In the scale-invariant case where Ug = 1, /Cn only 
depends on ki/k^ and /c2/^3, and lCn{ki, k2, k^) = 1 for all ki = k2 = k^. 

To allow for potentially divergent shapes, we allow the powers {p, r, s} to be negative as 
well as positive, and to make each /C„ nearly scale-invariant, we require the powers to satisfy 
p + r + s = 0. 

Each shape has a well-determined behavior in the squeezed limit, 

1 

JC7= Y.d^^x^^ + 0{x%). (2.19) 

m=R 

The set of {/C„} with the allowed combinations of {p, r, s} is given in Table 1 along with 
their divergence properties. The basis modes are also written in the equivalent short-hand 
notation used by [62]. Since we will use this notation to describe non-separable shapes in the 
next section we summarize it here: 

Kp = t: kf, = ^ E E ^I' Kp'^r = E fcf E E (2.20) 

i=l ^'^ i=l j^i i=l jy^i iy^ij 

with 

= 1 -|- dpq, ^pqr — ^pq{^qr ~l~ ^pr)- (2.21) 



2.4 Application of the basis to shapes arising in inflationary theory 

In this subsection, we illustrate the efficiency and accuracy allowed by our basis in describing 
shapes in the literature. First we discuss cases involving shapes and templates which are 
exactly expressed in terms of, or well-approximated by, linear combinations of the first 3 
modes of the basis. Then we extend the basis to include more divergent modes, and present 
the basis of shapes we will use under different divergence priors. 

2.4.1 Shapes exactly expressed in terms of {/Co — /C2} 

Some commonly considered templates can be exactly expressed in terms of the first three 
modes of the basis {/Cq ~ ^2}- 

The local shape, Siocal = can be derived from a simple ansatz for describing the 
nonlinear contribution to the primordial curvature perturbation in real space as a local effect 
[14], 

$(x) = <I>i(x) + /jvL ($i(x) - (<i>i(x)» . (2.22) 

Local shapes arise out of single-field slow-roll models, though the amplitude of the bispectrum 
in this case is predicted to be undetectably small [61, 63]. Large, local non-Gaussianity is 
predicted by a wide variety of other physically- motivated models, such as multifield inflation 
(e.g. curvaton scenario) [64, 65], (p) reheating mechanisms [66], and ekpyrotic inflation [65, 67]. 



-8- 



R = 
min{p, r, s) 




{p, r, s] 




K 


-3 


-2 


-1 








/Co 


(0,0,0) 








1 


Kqqo 


-1 


/Ci 


(-1,0,1) 
(-1,-1,2) 






1 
3 
2 
3 


1 

3 


K12 
Km 


-2 


/C3 
/C4 
/C5 

/Ce 


(-2,1,1) 
(-2,0,2) 
(-2,-1,3) 
(-2,-2,4) 




1 
3 
1 

3 
1 

3 
2 
3 


1 

3 


1 

3 


Ksa 

K222 
K24 
K222 

Kl5 


-3 


/C7 
/Cs 

/Cg 
/Cio 

/Cii 


(-3,1,2) 
(-3,0,3) 
(-3,-1,4) 
(-3,-2,5) 
(-3,-3,6) 


1 

3 
1 
3 
1 
3 
1 
3 
2 
3 


1 
3 


1 

3 


1 

3 


-* ^40 
K333 
K36 

K333 
K27 

""if 
K333 




etc. 















Table 1: The set of /C„ with the allowed combinations of {p, r, s} that satisfy p + r + s = 0. 
The ordering of the modes is according to their divergence behavior in the squeezed limit. 
The coefficients of their divergent terms in this limit, c?„m, are trivially related to the values of 
{p, r, s}, but are given here for convenience. We also give the equivalent short-hand notation 
used in [62] and summarized in (2.20). 



The constant shape, Sconst = 1 = /Co, was originally studied for its very simple form [62]. 
More recently the shape has been studied in the context of shapes arising from quasi-single 
field inflation (QSFI) models [23, 24, 41, 68]. The more general shape of QSFI is discussed in 
more detail below. 

Models with higher-derivative kinetic terms and/or non-trivial speeds of sound in the 
inflationary Lagrangian generally produce non-separable shapes, sensitive to the sum kt = 
ki + k2 + ^3 in the denominator, and thus cannot be exactly written in terms of a separable 
basis. The equilateral template [17], SgguU = — 2/Co-|-6/Ci— 3/C2, is widely used as a template to 
detect evidence of such shapes. Examples include generalized single-field models [7, 19, 69, 70], 
fe-inflation [10, 71, 72], ghost inflation [73], DBI inflation [74, 75], single-field non-slow roll 
and bimetric theories [21, 76]. For example, the single-field and DBI shapes. 



Sdbi 

S single 



1 

Km 



(2.23) 
(2.24) 



have a high cosine with the equilateral template. 

A general, effective field theory of inflation, in which the inflation perturbations are 
parameterized by a scalar field, vr, can contain many distinct operators which in themselves 
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give rise to specific bispectra [18], such as 



27^222 ,^ 

•^^^ = 1? — lis' (2.25) 
Kiiik'l 

_ AKni 2Kn 2Kf, 11 SKuh 



While each can typically be well-described by the equilateral template, a linear combination 
of these picking out the differences between them can yield a very different shape. This 
realization led to the generation of the 'orthogonal' template, Sorth = — 8/Co + 18/Ci — 9/C2. 

While inflation derived from a Bunch-Davies vacuum can be written in terms of a plane 
wave basis with positive k modes, excited states that are not in the Bunch Davies-vacuum can 
have initial states with both positive and negative k. Models motivated by non-trivial vacuum 
states can produce shapes with denominators containing ki+k2 — k^ (and its permutations), 
rather than kt [27-29]. Unlike the equilateral and local templates, these shapes peak in the 
flattened configuration, when k^ = ki -\- k2- While this shape again cannot be reconstructed 
perfectly using separable basis functions, an ansatz proposed as a proxy to this shape can be 
given by Senf = — 3/Co + 6/Ci — 3/C2 [28]. The shape has zero amplitude at ki = k2 = k^, 
making the conventional normalization at this configuration unsuitable for this template. 

2.4.2 Shapes well-approximated by {/Cq — ^2} 

Non-Gaussian templates to describe single-field theories are not limited to equilateral and or- 
thogonal shapes. Fast-roll single-field non-Gaussian models [20, 21] retain the scale-invariant 
spectra but relax the condition for slow-roll inflation. [56] showed these can be written in 
terms of seven constituents, four of which are Siocal, Sconst, ^1, and Sgingle- The remaining 
three constituent shapes are ^ 

^22 c- -^23 c, 
Knikt Kiukf Kmkf 

all of which have significant cosines with the local template. 

Other shapes exist in the literature that, while not separable, to some degree interpolate 
between the templates discussed above and hence can be reasonably-well described by linear 
combinations of {/Cq ~ ^2}- For example, non-Bunch-Davies vacua generate shapes that can 
be equilateral, local, or enfolded [29]. 

Quasi-single field (QFSI) models [23, 24, 41, 68] motivated by string theory and super- 
gravity inspired inflation contain multiple fields, but the extra fields have masses comparable 
to the Hubble scale. These models can be well described by a family of bispectrum templates 
dependent on a single parameter, u, 



3kik2k3Y^^ [8kik2k3k^ 
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where is the Neumann function of order i^. This shape interpolates between the constant 
and local templates. Another set of models that combine multiple fields and higher-derivative 
terms [30-32] also generate configurations that interpolate between standard shapes, spanning 
the local and equilateral templates. 

^Our shape, S„, is related to Battefeld and Grieb's A„, by Sn ~ An/Km. 
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1 2 1 2 



Figure 1: The application of the separable basis to describe shapes motivated by theoretically 
distinct models that span a wide array of configurations in A:-space. The local, equilateral, 
and orthogonal templates are explicitly separable, and are included only for reference. The 
other shapes are not separable but are approximated by templates using linear combinations 
of {/Co)^i)^2}- We construct an orthonormal basis, TZn, in /c-space with uniform weighting, 
using a Gram-Schmidt decomposition for /C^, for < n < 2, starting with n = . [Left panel] 
The cosines between each shape and the TZn- [Right panel] The cumulative cosine between a 
constructed template Y17=o ^i'^i ^^'^ the true shape. 

We use the basis modes, {/Cq — IC2} to create templates for these non-separable shapes, 
S's-S, Sdbi, Ssingie, and Sqfsi{'^)- To demonstrate this, we generate an orthonormal basis 
{TZn} using the Gram-Schmidt algorithm in fc-space, taking TZq = ICq, and create a template 
^template = Yli=o oinT^n as in (2.12) that reduces the covariance between the shape and tem- 
plate. The effectiveness of the template's fit can be quantified by the cumulative cosine. In 
Figure 1 we show how the shapes discussed above can be well modeled by templates using 
linear combinations of the {/Co — /C2} templates. In each case the cumulative cosine for the 
template and shape exceeds 0.98. 

2.4.3 Shapes well-approximated by more divergent {/Co — /C„} 

There are two strong motivations to extend template design beyond these three core templates. 
Firstly, expansions using the first three templates do not necessarily ensure that theoretical 
priors on the divergence properties are satisfied by the template. An example of this is the 
consistency relation that requires shapes of single-field inflation to vanish in the squeezed 
limit [61]. However, the orthogonal and enfolded templates constructed to describe single- 
field shapes tend toward a constant value in the squeezed limit. [18] proposed an orthogonal 
template, Sortho{2)j and [77] an enfolded template, Sf,nf{2)i that are somewhat more complex, 
using linear combinations of shapes that diverge as Xsq~'^, but they have the benefit of showing 
the correct divergence properties and more accurately reproducing the original non-separable 
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r 


-2 


-1 





1 


-1 




^2 = Siocal 


/Co 

2/Ci — /C2 


-2/Co + 6/Ci + 3/C2 = Seguii 


-2 




2/C5 - /Ce 


2/C4 — /Cg 


/Co + 3/C3 — 3/C4 
/C2 + 2/C3 — 2/C5 
2/C3 - /Ce 



Table 2: Shapes constructed from basis modes with maximal divergence [R < 0) which, 
through cancellations of the divergent terms, have a squeezed limit that diverges as x^g. These 
represent an irreducible set of component shapes, for each value of R, from which general, 
scale invariant, separable shapes can be constructed. 

shape. They can be written in terms of the /Cn modes as 

Sortho{2) = + P)Sequil 
Senf{2) = (1 +P)Sequil 

where p is a variable chosen to maximize the template's fit to the physical shape. 

Using our basis we can generalize this approach and write down classes of templates, 
denoted S^^^^t^^ constructed from basis modes with maximal divergence R that in the squeezed 
limit diverge as x^^, where r > R. In general, a shape written in terms of the basis will have 
a squeezed limit behavior given by 

1 

Ssq — O^n ^ ^ dnm^sq ~^ '-^(■^sq); (2.31) 
m=R 

with dnm summarized in Table 1. We find S'j/j can be written in terms of an irreducible set 
of shapes given in Table 2 for which andnm = for R < m < r. 

Siocal and Sequil are the only shapes constructed from R = —1 modes that respectively 
have —1 and vanishing divergence. There are an infinite set of shapes, however, with constant 
divergence described by /3/Co + (1 — /3)(2/Ci — IC2) where /3 is free parameter which could take 
any value except /3 = —2, for which the equilateral template is recovered. Instead of varying 
the parameter /3, we could instead select a value of f3 to generate a template from the set. 
/3 = —8 corresponds to the orthogonal template chosen by [18] to maximize the resulting 
shape's orthogonality with Siocal and ScquU- We could then choose to write general shapes in 
terms of linear combinations of {Sequil, Sortho, Siocal}, rather than OnlCn, 

^[-1,0] = OlESequil + OloSortho, (2.32) 
_i] = aESequil + OoSortho + OiLSlocal- (2.33) 

If these are the only shapes being used, the normalization constraint S'j/j _r] (A^O; ^Oi ^o) = 1 
fixes one a coefficient. 



P 



g^o + - 2/C2 + 



y /C4 + 3/C5 



'6^ 16^ 
P\ -/Co + — /C3 



-/C4 + -/Ce 



/Ce (2.29) 



(2.30) 
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We can extend this approach to include basis modes that diverge as 3y (J 



S'[_9 11 = OtESequil + Oto{Sortho + 6/C4 — 6/C3) + aiiSlocal + 2/C3 — 2/C5) 

+{l-aE-ao-aL){2IC3-ICe), (2.34) 

OtESequil + OoSortho + a.L{Slocal + 2/C3 - 2/C5) + /33(2/C3 - /Ce) 

+(l-/33-aL-a£-Qo)(2/C4-/C6) (2.35) 

OtESequil + OioSortho + OiliSlocal + 2/C3 — 2/C5) + /33(2/C3 — /Ce) + P^C^fCi — fCe) 

+(1 - /33 - /34 - - aij - ao)(2/C5 - /Ce), (2.36) 

OtESequil + CHoSortho + OtliSlocal + 2/C3 - 2/C5) + fizC^lCz " ^g) + /34(2/C4 - /Ce) 

+/35(2/C5 - /Ce) + (1 - /33 - /34 - /35 - - - ao)/C6. (2.37) 



Si 



-2,1] 
[-2,0] 
-2,-1] 
-2,-2] 



To tie this general approach to specific shapes in the literature, Sgj.tho{2) 'S'en/(2) can be 
written in this form by the following choice of coefficients: 



^ortho{2) — (1 +P)Sequil - pS[-2,l] 



'en/(2) = (1 +P)Sequil - pS[-2,l] 
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:,OiL 



, ao 



:,aL = 



(2.38) 



(2.39) 



The inclusion of extra basis shapes can be particularly important when the shape has 
undulations and is not just a smooth monotonic function. Shapes arising out of Galileon 
inflation are a good example of this. Imposing a Galilean symmetry on a single-field inflation 
model [77-80] gives rise to a non-Gaussian shape generated by three cubic interaction terms 
in the inflaton Lagrangian. While the shapes associated with each of these three operators, 
individually, are well-approximated by Sequil and Senf(2)i there exist combinations of them for 
which the resulting Galileon shape has little overlap with any of the shapes we have mentioned 
so far. Non-separable templates for Galileon inflation have been developed in [77] and [80] 
which have high cosines both with the underlying shape and each other. 

For illustrative purposes, we consider the shape presented in [77], based on equations 
(26)- (28) of this reference. When we use the Gram-Schmidt decomposition to construct a 
template with only the first three modes, we find a poor fit with a cumulative cosine of only 
0.13. The Galileon shape derives from a single-field action and a Bunch-Davies vacuum so 
theoretical consistency requires that it vanishes in the squeezed limit. IMotivated by this, if 
we fit the Galileon model using the 4 shapes in S'[„2,i]) we obtain a template with a cosine 
of 0.93. This reconstruction is not improved if we allow an unconstrained combination of the 
seven /Co — A^e modes. 

We can extend our approach to i? = —3 modes, and for example consider the following 
general shape that vanishes in the squeezed limit: 

S[--i^l\ = aESequil + Oto{Sortho + 6/C4 " 6/C3) -|- OL^Siocal " 2/C5 -|- 2/C3) -|- ^3(2/C3 - ICq) 
+/37(2/C7 - /Cii) + MSortho + 6/C8 - 6/C7) + MSiocal " 2ICg + 21Cj) 

+(1 -aE-ao-aL- P^- Pi- Pi- ^^fi " 2/Cio + 2/C7). (2.40) 

Fitting these eight distinct shapes in 'S'[_3 1] to the Galileon shape, we obtain an improved 
template with cosine of 0.99. 

The second reason to consider a basis including more divergent terms is that some 
inflationary scenarios, such as excited initial states and warm inflation, in which inflation 
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occurs in a warm radiation bath [25, 26] (see [81] for a review), can give rise to shapes that are 
more divergent than the local shape, with an overall divergence of Xsg~'^- This would suggest 
using an unconstrained combination of ICq — /Cg modes, or using constrained combinations of 
the R = —3 modes for which the xj^ divergent term vanishes. One such example of this is a 
template for warm inflation proposed by [56], 

Swarm = 1^2 + 1^7 — 1^9- (2-41) 

The realization that the differences between similar shapes can be important and provide 
an additional insight into the underlying model, implies that we should not just compare a 
small number of templates to the data. It is reasonable to extend beyond this and create 
more refined templates, sensitive to more than just properties that models have in common 
with the equilateral, orthogonal, and local templates. 

2.5 Shape classification and depiction 

The models discussed in the previous section reflect only a sample of the wide range of 
non-Gaussian inflationary shapes arising in the literature. Putting a coarse filter on their 
properties, one might characterize them using three descriptors: i) their divergence in the 
squeezed limit, ii) how many modes it takes to accurately describe them, and iii) the "family" 
to which they belong. 

Many of the physical shapes tend to be grouped in terms of a "family" resemblance 
to an existing template, reflecting the type of configurations of triangles with sides lengths 
ki,k2, and ks where the shapes have most of their power [62, 82]. For scale invariant shapes 
this is equivalent to studying the distribution of power over the space {f^, ^} for a fixed 
k^ > ki,k2- This space can be pictorially represented by a triangle with sides < ^ < 1 and 

1/2 < || < 1. We introduce it here in the context of the shapes already discussed, because 
we use this format to present some of our forecasting results. 

In Figure 2 we show examples of the shapes discussed in the previous section. Sconst = K-o 
is the archetypal component of a family with similar power over all scales, homogeneous over 
the whole triangular region plotted. "Squeezed" shapes have a bispectrum amplitude that is 
peaked in the top left-hand corner of the plot where ^ 1 and ^ = 1, while "equilateral" 
type shapes peak in the top right-hand corner where ^ = ^ = 1. "Flattened" shapes peak 
along the left edge, where ^ + ^ = 1. 

Of the shapes we've discussed so far, some clearly fall within these family categories: 
Siocal, Swarm, 5*4 and ^5 are "squeezed" shapes, while Sequil, Sdbi, Ssingle are "equilateral' 
and Sf>nf is "flattened". 

There exist other additional shapes generated by modes /C3 through /C7. For example. 
Figure 2 also includes three shapes that contribute to 5'[_2,i] that could describe a general 
single-field model with Bunch Davies vacuum. While each vanishes in the squeezed limit 
by construction, we find they differ from the equilateral shape in still having a component of 
their signal focused along the flattened configuration. The comparative size of this component 
correlates with the divergence of the shapes from which they are created, Sortho, Sio^aU and 
/C3. 

There are shapes that do not fall clearly into any of these families: Sortho{2) peaks in both 
the flattened and equilateral configurations, excited states can peak in squeezed and flattened 
configurations, and Sqsfi shapes interpolate between constant and local properties. Beyond 
this there are shapes with distinct undulating forms, the ScalUeon shape for example, that do 
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Equilateral 



{^) Slocal 
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(b) S, 



equil 



*1 

(C) Senf 



Quasi-single (y-1. 3) 




a.o 0.2 



ki/k2 

(d) Sortho(2) 




0.0 0.2 0.4 0.6 0.8 1.0 



(e) Sqsfi{u^I.S) 




0.4 0.6 
ki I (t.-i 



(f) SGalileon 




(g) Sortho + 6/C4 - 6/C3 




0.6 o.s 



(h) Slocal — 2IC5 + 2JC3 




Figure 2: Plots showing the comparative spatial distribution of non-Gaussian shape, 
5(^1, A;2, fcs), as a function of ki/k^ and k2/k^- From left to right we show [top] the lo- 
cal, equilateral, and enfolded separable templates, [middle] the orthogonal(2) template, and 
non-separable shapes derived from a QSFI model with v = 1.3 and Galileon inflation, and 
[bottom] shapes contributing to 'S'[_2^i] • All shapes are normalized to unity at the equilateral 
configuration = || = 1). The color scales for all but the local and QSFI shapes are the 
same to aid comparison. 
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not peak at either edges or corners. Moreover, not all shapes within each family are alike. 
For example, the local and warm shapes both peak in squeezed configurations, but their 
divergence properties in this region are different, leading to a low cosine between them. 

Given the breadth of bispectrum shapes that could be created, and the comparatively 
loose characteristics on which "families" are formed, there is strong motivation to ask how 
much information we can discern observationally about bispectra. This will help quantify how 
well we might determine the underlying non-Gaussian shape, if a detection of non-Gaussianity 
is made. 



3 Forecasting constraints on general shapes 

In the following analysis, we apply the separable, divergent basis and template classes from 
the previous section to assess how we can constrain the shape of primordial bispectra with 
upcoming CMB data. Our goal is to quantify what properties of shapes are measurable, and 
the respective roles of the experimental uncertainties and theoretical priors on determining 
distinguishability. 

Motivated by a broad cross-section of models in the literature, we will focus on shapes 
described by basis functions {/Cq —fCe} that are nearly scale-invariant and contain terms that 
are potentially as divergent as in the squeezed limit. 

We describe the Fisher matrix approach we use assuming a Planck-like CMB experiment 
in section 3.1. In section 3.2, we present the results of a principal component analysis for the 
set of shapes S[_2^r] with different divergence criteria in the squeezed limit imposed. Doing so 
generates the experiment's preferred orthogonal basis of best to worst measured bispectrum 
configurations, the principal components (PCs) and their corresponding uncertainties, subject 
to the theoretical prior. We consider the implications for shape normalization and the best- 
measured /c-configuration in sections 3.3 and 3.4, respectively, and finish in section 3.5 by 
quantifying our potential ability to determine and distinguish shapes. 

3.1 Fisher matrix approach 

We compute the 7x7 Fisher matrix for the amplitudes of the basis modes, fCn, {Inl'''^ ~ 
0, ...,6} defined in Eq. (2.16) as 

FifhLJh.) = fs.y E E ^^(^--^)?Sr'%^' (3-1) 

where {ahc} and {xyz} each sum over the 8 possible temperature (T) and polarization [E) 
combinations of bispectra: TTT, TTE, TET, ETT, TEE, ETE, EET, EEE. 

Given a general primordial shape expanded on the {/Cn} basis as in (2.16), the corre- 
sponding CMB reduced bispectrum is 

n 

where l^i^j^is ^^^^otes the reduced bispectrum of the basis function /C„ in (2.17)-(2.18). Here 
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we compute it as 



abc{n) 



1 



x^dx 



/,7(x)/^;(x)/,-(x) + {pr.}p 



An) ,2(n,-l) 

dkkP A1{k)je{kx 

TT ju . 



erms 



(3.3) 
(3.4) 



where p' is defined as in (2.18). 

We have modified the CAMB code [53] to numerical evaluate the values of /"^ and 

(n) 

then written code to appropriately combine them to form each /C^^^^^^. Specifically, we take 
kmin = 6.6 X 10~^ Mpc~\ kmax = 0.56 Mpc~^, and Xmax = 16.5 X 10^ Mpc. 

The covariance matrix we use from [83] is 



{Cov 



-\\abc,xyz 



-l\ax 



{C 

+ic 
+(c-^)?f 



{C 



— l\bz(/^—l\cx I ( /-i—l\bx ( /-^—l\cz 



with 



r<TT r'TE 

-l\ax _ I '^i 

r<TE AEE 



-1 



(3.5) 

(3.6) 
(3.7) 



Nl"" is the 



Here j^^^y is the overall fraction of the sky observed, and we assume /g^y = 0.8 
instrument noise for a correlation between observables a and x. We model CMB noise by 
considering the three lowest frequency bands of the Planck HFI instrument for temperature 
and E-mode polarization, as described in the Planck Bluebook [84]. We assume each frequency 
channel has Gaussian beam profile of width OpwuM and isotropic noise with error in X = T, 
of ax- The noise in each frequency channel c is then given by 



Era"' 



-1 



(3.8) 
(3.9) 



Our fiducial flat ACDM cosmology is described by the following parameters, which are 
consistent with the latest WMAP 9-year constraints [2]: Ofe/i^ = 0.02258, Vlji^ = 0.1109, 
A^(fco) = 2.43 X 10^^, Ug = 0.963, and r = 0.088. As has been done in other recent Fisher 
forecasts on non-Gaussianity parameters, such as [85], we consider the uncertainties on the 
non-Gaussian amplitudes independent of the uncertainties in the fundamental cosmological 
parameters that also affect the power spectrum, as these are comparatively small relative to 
the uncertainties from the bispectrum shape functions [86]. For this initial analysis, we neglect 
the effect of imperfect measurements of the lensing signal [87, 88], secondary anisotropics [14], 
and inhomogeneous sky coverage/noise on the constraints (e.g. [51, 89]). 



^ http://camb.info 
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3.2 Fisher matrix results 

A general bispectrum can be expanded in terms of either /C" or the component shapes, {Sx}, 
in S[ji^r]i given in (2.32)-(2.37), 

B^{h,k2,h){kMsr ^ ^^^^ ^ ^ m.lCnih, h) = Y: fLSxih, h, h) (3.10) 

n X 

While the Fisher matrix we used based on S^ji^^] automatically includes the additional priors 
to constrain the divergence properties, these could also be introduced into the /C" Fisher 
analysis by using Lagrange multipliers to systematically impose each divergence constraint. 
The latter makes no assumption apriori about what linear combinations of the shapes given 
in Table 2 should have their amplitudes constrained. While we use the shape expansion 
in our discussion below, we investigated both approaches and found they led to consistent 
conclusions. 

We use the Fisher matrix in terms of /C" to construct Fisher matrices for the component 
shapes in 5'[_2,r] for ^ = ~2, —1,0, 1. In Table 3 we give the £-space correlation coefficients 
based on (2.11), but here weighted by the data covariance between pairs of the component 
shapes, Sx and Sy, 

CorreiSx, Sy) = -jfi^- (3.11) 

yi'XX^YY 

This gives a measure of the similarity of the component shapes based on how they are mea- 
sured by the survey, integrated over all i combinations. 

We find the similarity between pairs of the four basis shapes in S'[_2,i], each of which 
vanishes in the squeezed limit, are primarily related to the divergence of the shapes from 
which they are derived. Sequil and Sortho + 6/C4 — 6/C3 are very similar to each other, while 
S'/oca/ — 2/C5-I-2/C3 and 2/C3 — /Ce have a reasonable degree of overlap. Interestingly the 2/C3 — /Cg 
shape also has significant similarities with the shapes that diverge as xj^. This is derived 
from its strong signal along the configurations between squeezed and flattened configurations, 
as discussed in section 2.5. The shape with xj^ divergence constructed from the R = —2 
modes, 2/C5 — /C3, is highly degenerate with the local template; essentially this implies the 
two are indistinguishable from one another using the CMB data. 

The unmarginalized errors, cr{f^^) = Fxx , give the uncertainty in the measurement 
of a specific template if the underlying theory is known to be wholly described by that 
template. The covariance matrices obtained from inverting the Fisher matrices give the 
uncertainties on the amplitudes of the component shapes, (T{fNL), marginalized over the 
freedom allowed by each model. We summarize the results in Table 4. 

The covariance matrix in each case can be diagonalized to obtain the orthonormal eigen- 
vectors, 

ii = '^CixSx, (3.12) 

X 

and associated eigenvalues, which give the variances a'^{bi) in the amplitudes of the eigen- 
vectors. These then provide a way to rank the best to worst measured bispectra. Given this 
orthonormal basis, any general bispectrum may be expanded as 

/atl-S* = ^^jCi. (3.13) 
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Table 3: Correlation coefficients between shapes that diverge as These shapes are 

components in the general template classes, /S'[_2,r]) for r < n. 



The principal components obtained by diagonalizing the covariance matrix are not immedi- 
ately 'shapes' in the way we considered so far. They have unit norm with respect to the 
component shape basis, I'^ixl'^ = 1) rather than being normalized at the equilateral con- 
figuration, Ylx'^ixSx{ko,ko,ko) = 1. 

If we restrict the shapes to those described by the first three modes, marginalization does 
not significantly alter the constraints from the unmarginalized errors, i.e. the three common 
templates are essentially the principal components (PC) of the covariance matrix, with the 
eigenvalues showing that the more divergent the shape, the better it is measured. 

As shown in Figure 3, when extended to shapes constructed of seven modes, the cor- 
respondence between the PC's and divergence remains. We find that, in general, divergence 
in the squeezed limit, followed by a second divergence measure, corresponding to the signal 
near the flattened configurations, can be used as coarse indicators of comparative constraining 
power with the CMB. For the general shape without any additional divergence constraints, 
the best measured PC is almost completely composed of the most divergent shape, ICq. The 
second best measured PC has dominant contributions from Siocal and 2/C5 — /Cg with which 
it is very degenerate. While one might expect the orthogonal shape to be dominant in the 
next PC, we find that 2/C3 — ICq, with which the orthogonal shape has overlap, but which 
vanishes in the squeezed limit, is also significant. Although the xjg divergence is cancelled 
in 2/C3 — ICq, it still has significant signal as the squeezed configuration is approached, which 
boosts how well it can be constrained relative to SequU- 

3.3 Drawbacks of normalization at the equilateral configuration 

As stated earlier, the PC's as they are originally generated, are not shapes in the usual sense 
because they are not bispectra normalized at ki = k2 = = kQ. They have a unit norm 
in terms of the basis shapes. With this normalization, as usual in PCA, their eigenvalues 
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Table 4: The uncertainties on the amphtudes of the component shapes, in the general tem- 
plate classes 'S'[-2,r] > that diverge as x^^ in the squeezed limit. We give both the unmarginahzed 
errors, assuming the underlying shape is exactly described by the component shape, and the 
marginalized errors if we allow the shape to be a general linear combinations of components 
consistent with the prior on the divergence properties. 
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Table 5: Properties of the principal components for each template class 5'[„2,r] in terms of 
their component shapes. The properties in the squeezed limit is determined by the value 
of r. The table provides uncertainties, for a unit norm eigenvector, cT(6j), and an effective 
(7{f]\fL{ei)), when the eigenvector is normalized consistently at the equilateral configuration. 

quantify which combinations of the basis shapes are best and worst constrained by data, and 
the eigenvectors can be combined to create general shapes. 

We can convert a{bi) to an effective cr{f]\fL{ei)), corresponding to the amplitude of 
each eigenvector shape normalized in the conventional way, cr{fNL{f^i)) = W{bi)(ii{kQ, /cq, A;o)|. 
Table 5, gives the values of a{bi) and (j{fNL{(^i))- We quote the results when both temperature 
and polarization data are included. We find that the exclusion of the E-mode polarization 
from the Fisher analysis does not noticeably change the shape of the principal components, 
but does increase the eigenvalues by about a factor of ~ 1.5 — 2 across all eigenvectors. 

Normalizing our PC's at the arbitrarily chosen equilateral configuration allows us to 
compare them to other shapes consistently at one point in /c-space. (y{fNL) does not in itself, 
however, quantify a shape's overall variance across all k. An analogous situation arises in 
quoting uncertainties on the power spectrum amplitude from two different surveys, say a 
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Figure 3: The relative importance of each component shape in the principal components 
of the general shape with xj^ divergence constructed from the first seven modes, S[_2-2] 
[upper panel] and the shape which vanishes in the squeezed limit, 'S'[_2,i] [lower panel]. The 
principal components are ordered from best measured, i = 1 to worst measured. The best 
measured modes in S[_2-2] are highly correlated with the most divergent shapes, /Cg and 
S local- Though the 2/C3 — /Cg mode vanishes in the squeezed limit, it still has significant 
weight as it approaches the squeezed configuration, which leads to it being a better measured 
vanishing mode than Sequil- 

large-scale CMB survey and a galaxy survey. Both surveys could quote uncertainties at a 
common arbitrary scale, say fco = 0.05h/Mpc, but while this uncertainty might represent the 
best measured scale for the galaxy survey, it would grossly overstate the minimum uncertainty 
in the CMB survey, which is best measured at a much larger scale. 

It is entirely possible for a well measured mode to have a significant part of its small 
variance located in the equilateral configuration, while a poorly measured mode could have its 
lowest variance in the equilateral configuration but be poorly measured over other regions of 
A;-space. Indeed we find this to be the case, given that the best measured shapes have signal 
peaked near the squeezed, rather than equilateral, configuration. This means that <j{fNL) is 
not a useful measure in itself to assess how well a shape can be measured. The overall spread 
in uncertainties from the best to worst eigenvector is much reduced when normalized at the 
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(g) 67 

Figure 4: Configurations of tlie principal components for S[_2-2] , a general shape that can be 
as divergent the squeezed limit. The plots show the amplitude of the eigenvectors 

for the best ei to worst ej measured modes as a function of versus The principal 
components are each normalized to be unity at the equilateral configuration. The large 
amplitude of 64 arises only because the eigenvector has a small amplitude in the equilateral 
configuration; when normalized here, the regions of large signal, towards the squeezed limit, 
are also boosted. 
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equilateral configuration and can in some cases produce a switch in the ordering of the modes 
for (T^/nl) relative to that of a(bi). This does not present an inconsistency in the analysis, 
but simply demonstrates the perils of considering a normalization at an arbitrary scale. 

Figure 4 shows the variety of profiles in the 2-dimensional ^) space shown in the 
triangle plots. Given that the power spectrum we consider is not perfectly scale invariant, 
there is some, small dependency of the bispectrum amplitude on the value of fea, described 
by p' in (2.18). The spatial profiles, however, in terms of ^ and || are fes-independent. 

The gradients in the PC configurations reflect the rough ordering from squeezed to 
flattened to equilateral as the modes span from best to worst. The complementarity of the 
eigenvectors, reflected by the different directions of gradients of the signals in the configuration 
space, has implications for the location of the best measured configuration, we discuss in 
section 3.4. 

3.4 Best measured A;-configurations 

In the analysis that follows, we avoid splitting up bispectra into shapes and amplitudes, 
normalized at an arbitrary configuration. Instead we consider the overall constraints on the 
bispectrum, B{ki,k2,k3), itself up to the constant normalization, given in (2.2), JnlS = 
kfk^klB{ki,k2,k3)/N. 

The eigenmodes and eigenvalues from the PCA provide a way to compute an error on a 
general /c-space bispectrum. 

We can calculate the posterior distribution of the uncertainties on JnlS given the data, 
with a theoretical prior given by the eigenvectors {cj}. 



~ n 

p{fNLS\D)= I lldhp{fNLS\b,)p{bi\D), (3.14) 

i=0 

n 

p{fNLS\bi) = 6{fNLS - ^ biei{ki, k2, fcs)), (3.15) 

1=0 

Under this assumption of Gaussian errors this gives the commonly used result, 

aHfNLS{ki,k2,k3)) = Y,a\bi)e^{ki,k2,k3f. (3.17) 

i 

This equation for computing the error can be applied to each set of PC's generated for each 
divergence scenario in the previous section. The errors in the (^) ^) configuration space can 
be plotted and the best measured /c-configuration, and the associated uncertainty, calculated 
for each scenario. 

(^{InlS) varies only very weakly across slices in k^; its functional form can be divided 
into a dependence on (f^, ^) and a weak dependence on k^, going as k'^^"'" ^\ for fixed ^ 
and For our choice of theoretical priors on the model, (T{fNLS) decreases with increasing 
k^. This is because the noise scales as the signal for the near scale-invariant theoretical prior 
we impose. An alternative prior would give very different dependencies on ^3. For example if 
we were to remove the theoretical prior all together and model the bispectrum amplitude as 
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Figure 5: As in Figure 4, but allowing the configurations of tlie principal components for 
5'[_2,i], a general shape that vanishes in the squeezed limit. 



bins in /c, the only constraints on the model come from the observational uncertainties, and 
the noise would diverge exponentially on small scales. 

The weak dependence implies that the uncertainties at one A:3 reasonably reflect 
the overall uncertainties if one were to marginalize over k^. Figure 6 shows the error on the 
/cs = 0.01 Mpc~^ slice for three different divergence cases. The location of minimum (t^/nlS) 
comes from the sum of the eigenmodes that is weighted by each mode's error, which arises 



- 24 - 



out of the complementarity of the degeneracy directions of the PC's. We find the location of 
the best measured configuration is consistent for the scenarios that diverge as x~q through to 
a constant in the squeezed limit for R = —2. This location is not situated in any one of the 
corners of the triangle plot associated with squeezed, equilateral, and flattened configurations. 
Instead it is somewhat centrally located adjacent to the flattened edge. The best-measured 
configurations are located at « 0.35 and ^ ~ 0.80, with minimum (T{f]sfLS) ~ 21. For the 
vanishing divergence prior, the best-measured configuration approaches the squeezed limit, 
as we have required the noise to scale as the shapes, which go to zero there. We also find 
that for shapes constructed from the local, equilateral, and orthogonal templates {R = —1), 
the best measured location spans a degeneracy direction also along the flattened edge, with 
minimum cr{fNLS) ~ 20. 

While the best measured region, in which the error is a minimum, is useful in the 
absence of knowledge about the theory, the signal-to-noise ratio, for a given underlying model 
can also help determine a survey's ability to distinguish between shapes. We consider this in 
the following subsection. 

3.5 Shape determination and distinguishability 

We now turn to discussing a central question of the paper: given a detection of non-Gaussianity 
using a specific template, what can be confidently inferred about the true underlying shape? 
We have already considered this from one perspective in section 3.2 by considering the uncer- 
tainties in ascribing a detection using a template to the template's shape itself. If we allow 
for the possibility that a detection using a specific template could be detecting the compo- 
nent of another shape allowed by the theoretical prior we are considering, then the errors on 
shape determination can increase significantly, especially for shapes that do not peak in the 
squeezed or flattened configurations. 

In this section we approach the question of shape distinguishability from a second di- 
rection, considering the range of possible general shapes, under a divergence prior, that could 
create the detected template signal and fit the bispectrum data within some confidence range. 
Such analyses have already been considered in the context of specific models, for example, 
how well we might disentangle a QSFI model (e.g. in (2.28)) from SgguU or Siocal ^ ^ function 
of ly [41, 68]. Here we extend this approach to a more general shape, and consider what impli- 
cations a detection with one of the common templates has for general models. For specificity 
we consider a subset of general shapes consistent with S'[_2,i], 

Sgen = [l - ax - ay - az)Sequii + oixSx + aySv + azSz, (3.18) 

where Sx,y,z can be {Sortho + - 6K3, Siocai - 2-K'5 + 2i^3, 2K3 - i^e}- This is investigating 
a general set of single-field inflation models from which Sortho{2) hi (2.38) and S'e„/(2) in (2.39) 
are drawn. 

How large must a template signal fj^^ be to be confident that the signal is not from 
a different, more general shape Sgeri^ We set this distinguishable detection threshold to be 
(T{fjfl), the error on fj^^ for the template, marginalized over f§i, the amplitude of the 
general shape. The marginalized constraint is computed by inverting the 2x2 Fisher matrix 
for ifJfL, f^i)- Thus we are comparing two shapes, where one is a template, and the other 
is a general shape, in which ax, ay, and az parametrize the deviation from SequU- 

In the simplest case, we allow only ax to be non-zero, such that Sgen is a linear com- 
bination of two shapes, Sequil and Sx, that varies with one parameter. In Figure 7, we 
show O'ifjfi) for the local, equilateral, and orthogonal templates when Sx takes different 
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(a) Diverges as 1/fc^ in the squeezed limit. (b) Vanishing divergence in the squeezed limit. 

tr(./NLS) 




(c) Diverges as in the squeezed limit. 

Figure 6: Example contour slices of cf{fNLS) under different divergence constraints, [top 
left] 5'[_2 -2], [top right] 5'[_2,i], and [bottom] S'f^i 



forms. The minimum value of for each template across all values of a recovers the 

unmarginalized errors of each. A detected value of /^^*' must be larger to produce a la 
detection of the equilateral shape, as opposed to a more general shape with a 7^ 0, while 
jZocaZ jjgygj- }jag ]-,g much larger than the unmarginalized cr(/^£"') to favor the local model 
over this general shape, because Sx and Siocal are weakly correlated. 

To illustrate the use of Figure 7, for example, a detection of fj^f^ = 40, while greater 
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than the unmarginahzed error of 19, would only be sufficient to rule out a false la detection 
of Sgen with Sx = Siocal — 2/C5 + 2/C3 for — 13 < Q < 2. On the other hand, if f^f^ is 
detected to be larger than 47, then Sgen of the 1-parameter form would be disfavored, as 
c^ifpfj^) is smaller than this over all values of a. Models with Sx = 2/C3 — /Cg are most easily 
differentiated from S^quil because they have the lowest correlation with SequU- 

An application of comparing constraints on ax from two distinct templates is to test 
whether a given model is consistent with or disfavored by the data. If two template measure- 
ments each individually remain consistent with two non-overlapping regions of a-space, then 
it would be clear that modeling the underlying shape with a alone is not able to produce a 
viable model. This would be true for dual measurements of {f^^^ = QO, f^j^° = 40} for 
Sx = Siocal — 2/C5 -I-2/C3, since they would imply non-overlapping ranges of a, —1.4 < a < 0.6 
versus a < —13.2 or a > 1.5 to each be consistent with the data. 

We can extend the same analysis to a comparison between templates and a 2-parameter 
general shape by allowing both ax and ay to vary simultaneously, while az is fixed to 
zero. For example, Senf{2) is a specific template for which this is true. In Figure 8 we show 
f7(/^™') and cr{f^j^) over different choices of the 2-dimensional space and find that there 
exist degeneracy directions that are not fully captured by the 1-dimensional projections in 
Figure 7. We find that remains close to the unmarginahzed value in this case as 

well. 

In the most general 3-parameter model, we can ask the question of whether there is any 
area of this space corresponding to a general model that vanishes in the squeezed limit, with a 
significant enough overlap with the local template to require that a potentially detected 
be much greater than the unmarginahzed value of 3.5. If this were the case, then it may be 
that a local template detection cannot definitively rule out a general shape that satisfies the 
single-field consistency relation. However, we find that nowhere in the parameter space does 
the cr(/^£) become greater than 4.1, showing that a detection of the local template above 
this threshold would effectively rule out a general shape, vanishing in the squeezed limit, 
subject to the assumption that it can be written in terms of our basis in 5[_2,i]. The same 
distinguishing power is not present for Siocal if we allow a weaker prior given by S[_2-i]- In 
this case the significant cosine between Siocal and 2/C5 — /Cg, means we may never be able 
to confidently attribute a detection with Siocal to be definitive evidence that the diverging 
signal is unambiguously Siocal- A long shot could be to additionally look at the correlation 
of the bispectrum signal with 2/C4 — /Cg which is mildly negatively correlated with Siocal and 
essentially uncorrelated with 2/C5 — ICq. 

This last point raises an interesting application of our study: to ask if there are distinct, 
new templates that we might use to learn about the origins of a detected non-Gaussian signal. 
In the context of models described by the ffist three modes, /Cq to /C2, the local, equilateral, 
and orthogonal templates are almost perfectly aligned with the principal components. If we 
extend the templates to include /C3 through /Cg, however, we find these no longer represent 
the PC's. For example, what might be the best way to extend the template pool to search for 
signatures of single-field inflation models with Bunch-Davies vacua? In the context of r = 1 
shapes, 2/C3 — /Cg is well-aligned with the best measured PC and is only mildly correlated with 
the existing templates which would make it a reasonable candidate to add as an additional 
template. We show the resulting constraints on general shapes in Figures 7 and 8. The 
figures show that this template probes regions of the allowed a-space which the equilateral 
and orthogonal templates do not constrain in the same way. Thus it may be possible to 
combine constraints from the common templates and motivated choices of a small number of 
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new templates, like 2/C3 — /Cg) to probe the underlying shape of non-Gaussianity. 
4 Conclusion 

At the heart of this work is the discussion about how uncertainties quoted on shape detection 
are inherently dependent on the underlying assumptions made about the shape. While a 
detection of non-Gaussianity with any template will be extraordinarily transformative in our 
field, its interpretation, in what it tells us about the underlying shape, has to be considered 
carefully in terms of our underlying theoretical prior we impose. 

We have presented an approach for quantifying how well upcoming CMB temperature 
and E'-mode polarization data can determine the shape of primordial non-Gaussianity under 
minimal assumptions. We proposed a set of polynomial divergent basis functions, {Kn}, that 
are well-tuned to describing many nearly scale-invariant, smoothly varying, but potentially 
divergent shapes discussed in the literature. We find we need only three to seven modes to 
generate matched templates to describe a wide range of physically motivated shapes. In this 
sense, the divergent basis is more efficient than the polynomial basis used in previous studies 
(e.g. [51]). Each /C n in our basis is generally divergent, but linear combinations of the ICn 
can be constructed to have cancellations in the squeezed limit, thus creating templates that 
are less divergent (e.g. equilateral shape). For example, S^quU and (2/C3 — KLq) both vanish 
in the squeezed limit, but still have a low correlation because the latter has more power near 
the flattened and squeezed configurations. 

Using the {/Cn} it is straightforward to form template classes, •S'jR^r]) that have specific, 
common divergence properties in the squeezed limit. Each class is constructed from an irre- 
ducible set of shapes, that while constructed out of a basis sets with maximum divergence 
(i? < 0), through cancellations of divergent terms, have squeezed limit x^^ {r > R). The 
choice of R controls how many basis modes are used to develop the templates, e.g. R = — 1 
includes /Cq through IC2, while R = —2 uses /Cq through ICq. As R becomes more negative it 
allows templates to be refined and shapes with a broader set of features across configuration 
space to be modeled. 

The classes allow templates to be developed with priors that are well-motivated by 
theories: S^ji^^ represents the class of all single-field models derived from a Bunch-Davies 
vacuum, S'[^^_i] in addition includes all multi-field models that diverge like the local shape 
in the squeezed limit, and S[ji _2] is the most general class of shapes which includes models 
from non-Bunch-Davies vacuum initial states. 

While the constituent shapes making up each class have the same divergence proper- 
ties in the squeezed limit, away from this limit they have power weighted differently in the 
configuration space. For example, we discuss a new shape, 2/C3 — /Ce, used in ^[/j i], that 
has the same squeezed limit behavior as the equilateral shape but has £-space cosines with 
the standard equilateral, orthogonal, and local templates of 0.07, 0.80 and -0.26 respectively. 
While the divergent terms cancel in the squeezed limit, 2/C3 — ICq has significant power just 
away from the squeezed and flattened configurations that differentiates it from the equilateral 
shape, and leads to it being most similar (though only mildly) with the orthogonal template. 

An added benefit of using the divergent basis and template classes to consider general 
shapes is that it ties together the methods we use to search for evidence of shapes with CMB 
data to LSS constraints from a scale-dependent halo bias, which probes the squeezed limits 
of shapes. It is well-known that templates for physical shapes which work for generating 
CMB predictions can fail when used for LSS predictions [42], because while CMB constraints 
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represent a weighted average over all fe-space configurations, the halo bias traces the squeezed 
limit region of /c-space only. Thus our approach provides a way of generating templates that 
can potentially be used consistently for both CMB and LSS studies. 

We adopt a Fisher matrix approach modeled on a Planck-like survey to estimate uncer- 
tainties on the amplitudes of shapes within each shape class, r. As summarized in Table 4, 
we computed the uncertainties on shape attribution under each prior and how these uncer- 
tainties on confidently being able to determine that a template is the true shape can change 
substantially dependent upon the type of prior we impose. 

We find that the best measured shapes are those with the strongest divergence and 
with principal power near squeezed and flattened fc-configurations. Though the conventional 
approach is to quote constraints at the equilateral configuration, ki = k2 = k^, we show, as 
summarized in Table 5, that this convention can mask how well or badly a shape is measured, 
as doing so has the effect of re-normalizing constraints such that badly measured modes can 
appear to have constraints similar to the best measured mode. 

Using the PCA results, we map out the fc-dependence of the constraints for a general 
shape given a prior, and show its dependence on the prior. For all but the r = 1 case, the 
best measured location is not in the equilateral configuration where shapes and constraints 
on /tvl are typically normalized, but in a configuration that is neither squeezed, flattened, or 
equilateral, but somewhere in between. This best measured location at roughly ki/ks w 0.35 
and k2/k^ ~ 0.80 arises out of the complementary gradients of the power in the PC's. For 
the r = 1 case, the best measured location is weighted more strongly towards the squeezed 
configuration, reflecting that the signal and the noise, with which it is correlated, both go to 
zero in this limit. 

Given our parametrization of a general shape under a divergence prior, we then ask 
how well it could be constrained using measurements of amplitudes of common templates, 
like the local, equilateral, and orthogonal templates. We focus on the class of general shapes 
that can represent the possible range of single-field models that vanish in the squeezed limit 
(r = 1). We calculate bounds on the subset of shapes that can remain consistent with 
constraints on the local, equilateral, and orthogonal templates, and find again-consistent with 
what we found earlier in the analysis-that templates with more power in the squeezed and 
flattened configurations provide more stringent constraints on this class of general shapes. 
Thus, the local, equilateral, and orthogonal templates serve different roles in constraining 
general shapes; the local template, if detected with sufficient amplitude, will rule out any 
shape of this type, while the equilateral and orthogonal templates serve to put constraints 
around different regions of the parameter space. In this sense, constraints from different 
templates can be complementary. 

Furthermore, a general (unknown) shape, will have different overlaps with the templates, 
creating a possibility that by combining constraints on templates, the overall constraint will 
shed more light on the underlying theory than any one constraint alone. We find it can also 
be advantageous to look for signals with a new, distinct template, beyond the three standard 
ones, that could help constrain models more efficiently; we explored the potential for using 
2/C3 — /Ce in this context. 

In this initial study we use somewhat idealized assumptions focusing on the effects of 
cosmic variance and Gaussian noise from a homogeneous sky coverage. We recognize the 
rich potential for further study to other basis sets, that better characterize sharp or os- 
cillatory features in bispectra, the presence of isocurvature modes, and stronger deviations 
from scale-invariance. To confidently attribute a primordial source to any measured non- 
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Gaussianity one would also want to fully account for contributions from astropliysical and 
instrumental sources, including gravitational lensing, inliomogeneous sky coverage, and sec- 
ondary anisotropies from astropliysical foregrounds. There is also the substantial question of 
how large-scale structure measurements, with sensitivity to the squeezed limit, can comple- 
ment the CMB data in constraining these general shapes, as well as whether 4-point statistics 
and checks of non-Gaussian consistency ansatzes can play a role. We are tackling some of 
these intriguing issues in work in preparation. 
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Figure 7: Detection thresholds on the amphtude of templates, fj^i, for distinguishing be- 
tween the template and a general shape, S = (1 — a)Sequal + oiSx^ at the la confidence level. 
[Top] Sx = Sortho + 6/C4 - 6/C3 [center] Sx = Siocai - 2/C5 + 2/C3, [bottom] Sx = 2/C3 - K.q. 
Blue, orange, red, and black curves denote f'j^i^\ Inl^^ InT^^ ^'n^l'^'' ^ respectively. 
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Figure 8: Detection thresholds on the amphtude of templates, [top] , [center] , 

and [bottom] f'pj^i~^'', for distinguishing between each template and two forms of a general 
shape S at the la confidence level. The general shapes considered are [left panels] Sgen = 
{I - ax - aY)Sequii + axiSortho + 6/C4 - 6/C3) + ariSiocai - 2/C5 + 2/C3) or [right panels] 
Sgen = {l-ax- aY)Sequii + oix{Siocai " 2/C5 + 2/C3) + ay (2/C3 - /Ce) . Contours for are 
not pictured, because the marginalized (t(/^£) remains close to its unmarginalized value over 
these 2-dimensional spaces. The case where Sx = Sortho + 6/C4 — 6/C3 and Sy = 2/C3 — /Cg is 
not pictured because Sx and 5y are nearly uncorrelated, thus no additional information is 
revealed beyond that in Figure 7. 
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